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Abstract— Repeated mutations in the SARS CoV-2 have resulted in the emergence of various life- 
threatening variants including omicron, triggering a systemic crisis in global human health. Due to the 
unavailability of any specific antiviral treatment, the effectiveness of the currently available drugs or 
vaccines against this variant remains doubtful. Phytochemicals from their known medicinal plants were 
selected with the antiviral, antifungal, anti-inflammatory, and antioxidant properties and screened to study 
their therapeutic potential and druglikeliness properties. A library of 15 phytochemicals from various 
medicinal plants was constructed and their binding energies were calculated by docking them against the 
spike protein of the SARS CoV-2 omicron variant using the AutoDock software. Subsequently, the 
Pharmacokinetic properties like ADME and drug likeliness properties were also calculated. 
Hydroxychloroquine (HCQ), an antimalarial drug was chosen as the standard in docking analysis, 
displaying the binding energy of -5.22 kcal/mol whereas, Glycyrrhetic acid showed the highest binding 
energy with the value of -9.02 kcal/mol. Besides, all these 15 phytochemicals showed higher binding 
energy (€ -7.00 kcal/mol) with the better pharmacokinetic properties leading them a viable drug candidate 
to treat the infection caused by omicron variant. 


Keywords— Drug likeliness, Insilico studies, Omicron, Phytochemicals, SARS CoV-2. 
I. INTRODUCTION Coronavirus (MERS-CoV) emerged in 2012, and the 


Severe Acute Respiratory Syndrome Coronavirus (SARS- 
CoV) in late 2003 [2]. 


The emergence and persistence of the coronavirus 
pandemic that was first identified in Wuhan, China in late 


December 2019, is still ongoing with a significant rise in 
the rate of mortality and morbidity. 


According to the World Health Organization (WHO) 
dashboard report, till 28th February 2022, 434,154,739 
confirmed cases have been reported globally, including 
5,944,342 deaths [1], and an increase in the number of 
cases has been observed day by day, thus raising the great 
public health concerns worldwide. 

COVID-19 is an infectious disease caused by SARS-CoV- 
2 (Severe Acute Respiratory Syndrome Coronavirus-2). 
Genomically, this virus is related to two known CoV- 
strains, namely the Middle East Respiratory Syndrome 
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The disease is characterized by symptoms such as fever, 
dry cough, headache, sore throat, fatigue, diarrhoea, and 
dyspnea followed by acute respiratory distress syndrome 
[3] and is predominantly known to spread from person-to— 
person [4]. In addition to this, the virus is also known to 
infect the upper and lower respiratory tract, liver, kidney, 
gut, heart, and nervous system and eventually lead to 
multi-organ damage [5]. Immunocompromised people 
with obesity, hypertension, diabetes, and cardiovascular 
disorders are found to be more susceptible[6]. 


SARS CoV-2 is an enveloped, positive-sense single- 
stranded RNA virus with a genome size of 26-32 kb [7][8], 
and is known to be transmitted between animals and 


196 


Kumar International Journal of Environment, Agriculture and Biotechnology, 7(5)-2022 


humans [9], belonging to the family Coronaviridae and 
genus Coronavirus [10]. 


The viral genome consists of 14 functional open reading 
frames (ORFs) including two noncoding regions at both 
ends followed by the multiple regions that encode for 29 
proteins including structural, nonstructural proteins (NSP), 
and accessory proteins. ORF la and ORF 1b encode 16 
NSPs (nspl-nspl6) such as Helicase (nsp13), RNA- 
dependent RNA polymerase (nsp12), Papain-like protease 
(nsp3), main protease (nsp5) also known as 3C-like 
protease (3CLpro), and 20 methyltransferase (nsp16), 
which are likely to be involved in transcription, 
replication, and pathogenesis and play a vital role in the 
life cycle of the pathogen. Whereas structural proteins 
include, Spike surface glycoprotein (S), a small envelope 
protein (E), membrane protein (M), and nucleocapsid 
protein (N), which are essential for viral assembly and the 
production of a structurally whole viral particle. And the 
remaining 9 are the accessory proteins that provide a 
selective advantage in the infected host cell [11] [12]. 


The S protein plays an integral role in viral attachment 
through the recognition of host cell receptors and promotes 
the fusion of the viral and host cell membranes, to 
facilitate the viral entry into the host cells. This protein is 
composed of two functional subunits: an amino (N)- 
terminal S1 subunit and a carboxyl (C)-terminal S2 subunit 
with a single transmembrane region anchor. The SI 
subunit consists of an N-terminal domain (NTD) and a 
receptor-binding domain (RBD), along with two C- 
terminal domains. Whereas the S2 subunit consists of a 
fusion peptide (FP), heptad repeat 1 (HR1), central helix 
(CH) region, connector domain (CD), heptad repeat 2 
(HR2), transmembrane (TM) region, and the cytoplasmic 
tail (CT). The cleavage site at the border between the S1 
and S2 subunits is called S1/S2 protease cleavage site [13]. 


In the native state, the S protein exists as an inactive 
precursor. During the infection, the S protein gets activated 
by the TM protease serine 2 (TMPRSS2), a type 2 TM 
serine protease located on the host cell surface, by cleaving 
the S protein into S1 and S2 subunits, that are required for 
activating the membrane fusion domain after the entry of 
the virus into the host cell [14]. Later, the S1 subunit is 
shed and the S2 subunit undergoes large conformational 
change compared to its prefusion state that mediates the 
membrane fusion process, thus, inserting the fusion 
peptide into the host cell membrane, followed by the 
release of viral RNA into the host cell [13], where the 
polyproteins are translated from the RNA genome. Protein 
cleavage and the assembly of the replicase-transcriptase 
complex, make the replication and transcription process of 
the viral RNA genome occur. The replicated viral RNA 
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and the structural proteins that are newly synthesized are 
assembled and packaged in the host cell before the release 
of viral particles [14]. 


Since, the S protein plays an active role in the recognition, 
attachment, and fusion of the viral particle into the host 
cell, thus designing an effective antiviral drug or inhibitor 
by targeting this protein can be the best solution to treat the 
infection. 


SARS CoV-2 virus, which was previously identified in 
Wuhan, China in late 2019, has undergone multiple 
significant mutations resulting in the new variants 
including Alpha, Beta, Gamma, Delta, Omicron, and so 
on. 


Out of these variants, Omicron is one of the most rapidly 
spreading SARS CoV-2 variants across the world and has 
been classified as a “variant of concern (VOC)" by the 
World Health Organization (WHO) [15]. 


The first case of the omicron variant (B.1.1.529.1), was 
reported in South Africa on 23 November 2021 [16] and as 
of 14 December 2021, it had spread to more than 76 
countries [17]. This variant is known to be the highly 
mutated strain compared to the other VOCs (Alpha, Beta, 
Gamma, and Delta) with 50 mutations accumulated 
throughout the genome. Out of which, at least 32 
mutations are found in the spike protein, which is twice as 
many as the Delta variant. These mutations affect the 
biological characteristics of the omicron variant, by 
increasing the transmissibility, causing immune escape, 
and enhancing the virulence. One of the Insilico studies 
showed that the infectivity of the Omicron variant might 
be more than 10-fold higher than that of the original virus 
(Wuhan) [18][19]. 


The emergence of this variant has posed an increased risk 
to the global public health, concerning the effectiveness of 
the currently available drugs or vaccines. At present, most 
clinical studies are in progress, to study the effect of 
antiviral drugs, corticosteroids, antimalarials, and 
monoclonal antibodies against this variant. 


Currently, only a few repurposed FDA-approved drugs are 
being used to treat the infection including Remdesivir, 
Lopinavir/Ritonavir, Oseltamivir, Favipiravir, Rifampicin, 
Letermovir, Azithromycin, Cholorquine, and 
Hydroxychloroquine [20]. Corticosteroids including 
dexamethasone [21]and methylprednisolone [22], are also 
being used to treat the infection. 


Along with this, scientists worldwide are studying the 
effectiveness of these drugs alone or with a combination of 
other antiviral drugs. However, few of the drugs showed 
certain side effects like cardiotoxicity, hematologic 
toxicity, hepatotoxicity, and nephrotoxicity. 
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E.g., Chloroquine and Hydroxychloroquine which were 
found promising earlier have resulted in neuromyopathy, 
cardiomyopathy, and retinopathy and was found to 
decrease in-hospital survival and increased the frequency 
of ventricular arrhythmias in hospitalized patients [23]. 


Serum neutralization assay against the omicron variant 
showed that the neutralization activity of the 
Pfizer/BNT162b2 mRNA vaccine declined the serum 
activity by 41.4 fold compared to the original strain 
(Wuhan) [18][24]. In addition to this, an in vitro study of 
imdevimab and casirivimab was found to be resistant to 
authentic omicron variant, which was effectively used to 
prevent Delta infection [25]. Thus, with the decrease in the 
vaccine efficacy and unavailability of any specific antiviral 
drug candidate against the omicron variant, the world is 
still battling to overcome the Pandemic. 


India is home to more than 8000 species of medicinal 
plants [26]. These plants have been known for their 
immense potential and tremendous properties against 
various infectious health-related 
complications from ancient times. However, these plants 
are a rich source of chemically active compounds 
including alkaloids, flavonoids, terpenoids, polyphenols, 
tannins, saponins, etc. possessing a variety of biological 
deliberations, including antioxidant, antibacterial, 


diseases and 


antifungal, anti-inflammatory, and antiviral properties 
[27]. 


At present, plant-based phytochemicals are attracting the 
focus of modern world healthcare researches, for 
developing a vaccine or treatment against the infection, 
since they are known to be less toxic with minimal side 
effects. It is reported that approximately 70-80% of 
mainstream medicines originated from natural products 
[28]. In addition to these, plant-based phytochemicals offer 
attractive, effective, and holistic drug action against 
pathogens. Thus, extracting the plant-based compounds 
and herbal treatments from these plants can be a cure and a 
solution against omicron infection. 


In this regard, Insilico studies were carried out to identify 
the effect of the phytochemicals that can play a role as an 
inhibitor by studying their binding affinity towards the 
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spike protein. In the modern world, Insilico studies are 
gaining much attention worldwide for their advanced 
strategies and effective techniques for identifying putative 
drug candidates that can be further led to in vitro and in 
vivo assessments. These methods basically provide the 
most accurate results, require shorter time duration, are 
cost-effective, and fall in the dry lab category, thus being 
widely used in understanding and predicting the drug 
ability in the early stages of drug discovery [29]. 


The current study is aimed to explore the drugs with 
a therapeutic potential action of the phytochemicals as an 
effective inhibitor against the spike protein of omicron 
variant, using Computational pharmacology and virtual 
screening techniques. A library of 15 phytochemicals was 
prepared from various medicinal plants, based on their 
therapeutic potential against infectious diseases. The 
binding energy of all these phytochemicals was calculated 
and compared with that of HCQ (an antimalarial drug) 
which was chosen as the standard reference drug for 
comparison. Additionally, ADME analysis and drug 
likeliness studies were also conducted for understanding 
the safety and efficacy of the phytochemical, which can 
play a role as the drug candidate for the further 
development of an effective drug or inhibitor against the 
infection caused by the omicron variant. 


II. MATERIALS & METHODOLOGY 
2.1 Preparation of receptor: 


The cryo-electron microscopic 3-Dimensional structure of 
SARS CoV-2 Omicron variant spike protein was 
downloaded from RCSB-Protein Data Bank with PDB ID: 
7T9J [30] and was visualized by using Chimera software 
1.14 [31]. In this study, Spike glycoprotein was chosen as 
the target site. Since, it plays a crucial role in its 
attachment, fusion, and viral entry into the host cell. The 
spike protein of the omicron variant consists of 3 chains, 
Chain A, B, and C as shown in Fig. 1. Out of which, Chain 
C (represented in Fig. 2) was selected for the docking 
process, eventually energy minimization and protein 
optimization was carried out by using the Swiss PDB 
Viewer [32]. 
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Fig.1: 3D structure of SARS CoV-2 omicron variant spike protein (PDB ID: 7T9J), Chain A (orange-red colored), Chain B 
(yellow colored), & Chain C (purple colored) 


Fig.2: 3D structure of Chain C of SARS CoV-2 omicron variant spike protein 


2.2 Protein optimization and energy minimization: 


This step is carried out to remove the additional molecules 
and the extra spaces from the pdb protein files, by 
arranging the protein coordinates into the proper order. 
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The pdb structure of the protein downloaded from the 
RCSB-PDB database was opened in WordPad and the 
individual chain (Chain C) of the protein was retrieved 
(represented in Fig. 3) and saved in ‘all files’ format. 
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Fig.3: Protein Optimization using WordPad 


Energy minimization of the protein chain was carried out 
by using the Swiss PDB Viewer tool (v4.1.0), followed by 
removing the additional swiss coordinates that were added, 
by using WordPad. 

2.3 Ligand preparation: 

A library of 15 phytochemicals from numerous medicinal 
plants was constructed, based on their therapeutic 
properties. The 3-Dimensional structures of all these 


phytochemicals including HCQ (PubChem Id: 3652), an 
antimalarial drug chosen as a standard reference drug were 


downloaded from the PubChem database in structure-data 
file (sdf) format. Later, the structure-data file (.sdf) format 
of these phytochemicals was converted into the pdb format 
by using Open Babel software 2.4.1 [33] and was 
visualized in Chimera software 1.14. The pictorial 
representation of the 3D structure of Glycyrrhetic acid 
visualized by using Chimera software is represented in 
figure 4. 


Fig.4: 3D structure of Glycyrrhetic acid. 


And Table 1 consists of the name list of all the 15 selected phytochemicals with their PubChem ID’s and the scientific name 
of the medicinal plant from which it is obtained followed by their therapeutic properties. 
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Table 1: List of phytochemicals isolated from few selected medicinal plants. 


Botanical name of the | Phytochemical PubChem Id Therapeutic properties 

plant 

Elodea canadensis B- Sitosterol 222284 Antioxidant [34][35], anti- inflammatory drug [36] 
and antiviral [37]. 

Syzygium claviflorum Betulinic acid 64971 Anti-HIV [38], antimalarial [39], anti-inflammatory, 
antioxidant and anticancer [40]. 

Desmodium gangeticum | Desmodin 13338925 Antioxidant [41] and used in the treatment of 
Alzheimer’s disease [42]. 

Momordica foetida Foetidin 15945065 Anti-diabetic [43]. 

Glycyrrhiza glabra Glycyrrhetic 10114 Antiviral [44], immunomodulatory, anti- 

acid inflammatory and hepatoprotective [45]. 

Ferula assa-foetida Kamolonol 46883037 Antibacterial [46] and cardioprotective property 
[47]. 

Marrubium globosum Marrubiin 73401 Antioxidant, anti- diabetic antinociceptive, 
antigenotoxic, antioedematogenic, antispasmodic, 
analgesic, gastroprotective, immunomodulating, 
cardioprotective and vasorelaxant properties [48]. 

Salvia tomentosa Maslinic acid 73659 Anti-inflammatory [49], anticancer [50], anti 
diabetic [51], antioxidative, cardioprotective [52]and 
neuroprotective properties [53]. 

Ophiopogon japonicas Oleanolic acid 10494 Antioxidative, antiviral, anti-inflammatory, 
immunomodulatory and cardioprotective properties 
[45][54]. 

Prunella vulgaris Oleanane 9548717 Anticancer, antimicrobial, anti- inflammatory, 
cytotoxic and hepatoprotective [55]. 

Gladiolus italicus Ursolic acid 64945 Antiviral [56], antioxidant, anti-inflammatory and 
anticancer [45]. 

Withania somnifera Withaferin A 265237 Anti platelet, anti herpetic, anti-inflammatory, 
anticancer, antileishmanial and immunosuppressive 
properties [57]. 

Withania somnifera Withanone 21679027 Antiviral [58], anticancer [59] and neuroprotective 
[60]. 

Withania somnifera Withanolide 53477765 Antiviral [61], anti-inflammatory 
Immunomodulatory and anti-cancer properties [62]. 

Plumbago zeylanica Zeylanone 5276618 Antimicrobial and cytotoxic properties [63] 


2.4 Prediction of Active Binding site 


An extensive literature survey was conducted to predict the 
active binding sites on the receptor molecule for docking. 
Few of the reports suggested that the amino acid residues 
such as R 493, S 496, and R 498 are considered the active 
binding amino residues present on the spike protein of the 
omicron variant [30][64][65]. 


2.5 Molecular Docking of Receptor and Ligand 
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effects of each of the 
phytochemicals were evaluated and compared with the std. 
HCQ by the docking method using AutoDock software 
1.5.6 [66]. Prior to the docking process, the protein and the 


The probable inhibitory 


ligands were prepared with the addition of polar hydrogen 
molecule and Gasteiger charges by removing the water 
molecules and a few torsion adjustments were made 
specifically for ligands and were saved in pdbqt format. 
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The dimension of the grid box was selected as 38À X 44À 
X 56À with the grid values of X - coordinate = 203.458, Y 
— coordinate = 181.765 and Z — coordinate = 269.208. The 
size and dimension of the grid box were selected based on 
the position of active binding sites of the omicron variant, 
with a total of 100 genetic runs, while the other parameters 
were set as default. During the docking process, the 
receptor (spike protein) was kept rigid and the ligand 
(phytochemical) was flexible. 


The binding affinities of these phytochemicals towards the 
targeted protein were identified and compared with HCQ. 
The data generated from the docking process showed how 
well the ligand gets interacted with the protein of interest. 
And these values were further analyzed by using MGL 
tools 1.5.6 [66]. 


Later, the output files obtained from the docking process 
were used to study the molecular interactions between the 
protein-ligand complexes including hydrogen bonds by 
using Ligplot+ 2.2 software [67]. 


2.6 ADME analysis & Drug likeliness 


The sdf structure of the phytochemical compounds 
downloaded from the PubChem database was fed into the 
SwissADME server [68] and converted into canonical 
SMILES format. These SMILES structures were used to 
study the ADME (Absorption, Distribution, Metabolism, 
and Excretion) properties of the screened phytochemicals, 
and the data obtained from this showed that how well a 
chemical drug (phytochemical) is processed in a living 
organism. 


Some of the important parameters related to the drug 
likeness properties of the compounds including molecular 
weight, number of hydrogen bond donors, number of 
hydrogen bond acceptors, and molar refractivity were 
computed utilizing Lipinski's rule of five [69]. 
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III. RESULTS 


The current study was focused on identifying the specific 
antiviral drug or inhibitor by targeting the spike 
glycoprotein of the omicron variant by drug repurposing 
technique since this protein plays a prime role in its 
attachment and internalization into the host cells. Thus, 
identifying a new drug or inhibitor which can inhibit the 
binding of the spike glycoprotein to the host cell can be a 
solution against the infection. 


3.1 Virtual Screening and Visualization 


Molecular docking is the most promising tool in the drug 
discovery process, which is used to predict the binding 
affinity of the ligand towards the targeted protein molecule 
when the 3D structure of the target protein is known. It is 
also used to study the interactions between the protein- 
ligand complexes at the atomic level and characterize the 
behavior of ligands in the binding site of target proteins as 
well as elucidate fundamental biochemical processes [70]. 
Because of their wide applications, this method is mostly 
preferred in all pharmaceutical industries. Hence, this 
method was employed to study the therapeutic potential of 
the selected phytochemicals including HCQ against the 
SARS CoV-2 omicron variant based on their binding 
affinity of the phytochemical (ligand) towards the targeted 
protein (spike protein). 


Higher the negative value of the binding energy stronger 
the binding of the phytochemical towards the targeted 
protein. 

The binding energy of all the 15 phytochemicals against 
the spike protein of the SARS CoV-2 omicron variant was 
calculated by AutoDock software 1.5.6. The 2D structure 
of the selected phytochemicals along with their binding 
energies and interacting residues are represented in Table 
2. 


202 


Kumar 


International Journal of Environment, Agriculture and Biotechnology, 7(5)-2022 


Table 2: The binding energy of the compound towards the spike protein of the omicron variant along with their interacting 
residues with the grid dimension of 38 X 44 X 56 (A) and grid value of X = 203.458, Y = 181.765 and Z = 269.208 ( The 2D 
structure of all the ligands represented in the above table are downloaded from PubChem database) 


Name of the ligand 


2D-Structure 


Binding 
energy 
(kcal/mol) 


Interacting amino acid residues 


Hydroxychloroquine 


e^ 


H 
o 


J 


E 


-522 


Hydrogen bond (H- bond): ARG 
403, ARG 493, ARG 498 and 
HIS 505. 


Hydrophobic interactions (H-I): 
TYR 449, TYR 453, SER 494, 
TYR 495, SER 496, PHE 497 
and TYR 501. 


B- Sitosterol 


-7.45 


H- bond: HIS 505. 


H-I: ARG 403, TYR 453, ARG 
493, SER 494, TYR 495, SER 
496, PHE 497 and TYR 501. 


Betulinic acid 


-7.46 


H- bond: ARG 493. 


H-I: ARG 403, TYR 453, SER 
494, TYR 495, SER 496, PHE 
497 and HIS 505. 


Desmodin 


-7.07 


H- bond: SER 494. 


H-I: ARG 403, TYR 449, TYR 
453, ARG 493, TYR 495, SER 
496, PHE 497 and HIS 505. 
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Foetidin 


-7.36 


H- bond: ARG 493. 


H-I: ARG 403, TYR 449, TYR 
453, SER 494, TYR 495, SER 
496, PHE 497, TYR 501 and HIS 
505. 


Glycyrrhetic 


acid 


-9.02 


H- bond: TYR 453 and SER 494. 


H-I: ARG 403, ARG 493, TYR 
495, SER 496, TYR 501 and HIS 
505. 


Kamolonol eo 


-7.41 


H- bond: TYR 453 and SER 494. 


H-I: ARG 403, ARG 493, TYR 
495, SER 496, PHE 497, TYR 
501 and HIS 505. 


Marrubiin 


-7.09 


H- bond: SER 496. 


H-I: ARG 403, TYR 453, SER 
494, TYR 495, PHE 497, TYR 
501 and HIS 505. 


Maslinic acid 


ISSN: 2456-1878 (Int. J. Environ. Agric. Biotech.) 
https://dx.doi.org/10.22161/ijeab. 75.20 


-8.17 


H- bond: ARG 493. 


H-I: ARG 403, TYR 453, SER 
494, TYR 495, SER 496, PHE 
497, TYR 501 and HIS 505. 
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Oleanolic acid 


-8.68 


H- bond: ARG 493. 


H-I: ARG 403, TYR 453, SER 
494, TYR 495, SER 496, TYR 
501 and HIS 505. 


Oleanane 


-7.63 


H-I: ARG 403, ARG 493, SER 
494, TYR 495, SER 496, PHE 
497, TYR 501 and HIS 505. 


Ursolic acid 


-8.54 


H- bond: ARG 493. 


H-I: ARG 403, TYR 453, SER 
494, TYR 495, SER 496, PHE 
497, TYR 501 and HIS 505. 


Withaferin A 


-7.53 


H- bond: TYR 449 and SER 496. 


H-I: ARG 403, TYR 453, ARG 
493, SER 494, TYR 495, PHE 
497, TYR 501 and HIS 505. 


Withanone 


-8.12 


H- bond: ARG 403, SER 494 and 
SER 496. 


H-I: TYR 449, TYR 453, ARG 
493, TYR 495, PHE 497 and HIS 
505. 
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Withanolide -7.13 H- bond: SER 494 and SER 496. 
H-I: ARG 403, TYR 453, ARG 
493, TYR 495, PHE 497, TYR 
501 and HIS 505. 

Zeylanone -7.21 H- bond: ARG 493 and SER 496. 


H-I: TYR 453, SER 494, TYR 
495, PHE 497, TYR 501 and HIS 
505. 


From Table 2, we can observe that the binding energy of 
all the 15 selected phytochemicals was higher (< -7.00 
kcal/mol) than the HCQ which was chosen as the std. 
reference drug to study and compare the binding affinity of 
these phytochemicals against the spike protein. 


The binding energy of standard drug HCQ was found to be 
-5.22 kcal/mol whereas, Glycyrrhetic acid showed the 
highest binding energy among all 15 phytochemicals with 
the value of -9.02 kcal/mol towards the active site of the 
target protein, followed by Oleanolic acid with the value of 
-8.68 kcal/mol, Ursolic acid -8.54, Maslinic acid -8.17, 
Withanone -8.12, Withanolide -7.73, Oleanane -7.63, 
Withaferin A -7.53, Betulinic acid -7.46, B- sitosterol - 
7.45, Kamolonol -7.41, Foetidin -7.36, Zeylanone -7.21, 
Marrubiin -7.09 and Desmodin which showed the lowest 
binding energy of value -7.07 kcal/mol. The binding of the 
receptor protein  (7T9J C) with each ligand 
(phytochemical) is represented in Figure 5(a)-20(a). 
Protein-ligand interactions of each phytochemical were 
estimated by Ligplot+ 2.2 software and are represented in 
Figure 5(b)-20(b). 


ARG 403, ARG 493, ARG 498 and HIS 505 were the four 
amino acid residues involved in the hydrogen bond with 
HCQ. TYR 449, TYR 453, SER 494, TYR 495, SER 496, 
PHE 497 and TYR 501 were the amino acids involved in 
hydrophobic interactions with the spike protein. 
Glycyrrhetic acid formed hydrogen bonds with the TYR 
453 and SER 494. ARG 403, ARG 493, TYR 495, SER 
496, TYR 501 and HIS 505 were the amino acid residues 
of spike protein that interacted with Glycyrrhetic acid by 
hydrophobic interaction. Similarly, Oleanolic acid showed 
the hydrogen interaction with the ARG 493 whereas, ARG 
403, TYR 453, SER 494, TYR 495, SER 496, TYR 501 
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and HIS 505 showed the hydrophobic interaction with the 
amino acid residues of the spike protein. ARG 493 amino 
acid residue of spike protein was involved in hydrogen 
bonding with the Ursolic acid. Besides, ARG 403, TYR 
453, SER 494, TYR 495, SER 496, PHE 497, TYR 501 
and HIS 505 were involved in hydrophobic interaction 
with the targeted site of the protein. Maslinic acid formed a 
hydrogen bond with the ARG 493. In addition to this, 
ARG 403, TYR 453, SER 494, TYR 495, SER 496, PHE 
497, TYR 501 and HIS 505 were the amino acids involved 
in hydrophobic interactions with the spike protein. 
Withanone showed hydrogen bonding with ARG 403, SER 
494 and SER 496 amino acid residues and moreover 
interacted with TYR 449, TYR 453, ARG 493, TYR 495, 
PHE 497 and HIS 505 amino acids residues by 
hydrophobic interactions. Withanolide formed hydrogen 
bonds with the SER 494 and SER 496 residues. ARG 403, 
TYR 453, ARG 493, TYR 495, PHE 497, TYR 501 and 
HIS 505 amino acid residues of spike protein interacted 
with Withanolide by hydrophobic interaction. Oleanane 
showed the hydrophobic interaction with the ARG 403, 
ARG 493, SER 494, TYR 495, SER 496, PHE 497, TYR 
501 and HIS 505 amino acid residues of the spike protein, 
whereas it didn’t show hydrogen interaction with the any 
of the amino acid residues of the spike protein. Withaferin 
A showed hydrogen bonding with TYR 449 and SER 496. 
Moreover, it interacted with ARG 403, TYR 453, ARG 
493, SER 494, TYR 495, PHE 497, TYR 501 and HIS 505 
amino acid residues via hydrophobic interaction. The 
Betulinic acid formed a hydrogen bond with the ARG 493. 
ARG 403, TYR 453, SER 494, TYR 495, SER 496, PHE 
497 and HIS 505 showed the hydrophobic interaction with 
the amino acid residues of the spike protein. 
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Similarly, B-sitosterol showed the H- bond with HIS 505, 
whereas, hydrophobic interactions with the ARG 403, 
TYR 453, ARG 493, SER 494, TYR 495, SER 496, PHE 
497 and TYR 501 amino acid residues with the targeted 
protein site. TYR 453 and SER 494 residues of spike 
protein were involved in H-bond formation with 
Kamolonol. Besides, ARG 403, ARG 493, TYR 495, SER 
496, PHE 497, TYR 501 and HIS 505 residues showed the 
hydrophobic interactions with the spike protein. Foetidin 
formed H-bond with ARG 493 amino acid residue. 
Moreover, Foetidin interacted with the ARG 403, TYR 
449, TYR 453, SER 494, TYR 495, SER 496, PHE 497, 
TYR 501 and HIS 505 residues hydrophobically. 


Fig. 6(a): 7T9J C - p- Sitosterol 
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Arg 403(C) 


Zeylanone showed hydrogen bonding with ARG 493 and 
SER 496. In addition to this, Zeylanone interacted with 
TYR 453, SER 494, TYR 495, PHE 497, TYR 501 and 
HIS 505 amino acid residues via hydrophobic interaction. 
Similarly, Marrubiin formed a hydrogen bond with the 
SER 496 amino acid residue. Furthermore, it interacted 
with the ARG 403, TYR 453, SER 494, TYR 495, PHE 
497, TYR 501 and HIS 505 residues hydrophobically. The 
Desmodin formed a hydrogen bond with the SER 494 
amino acid residue. In addition to this, Desmodin 
interacted with the ARG 403, TYR 449, TYR 453, ARG 
493, TYR 495, SER 496, PHE 497 and HIS 505 residues 
hydrophobically. 
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Fig. 15(a): 7T9J C — Oleanane 
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Fig. 13(b): 7T9J_C - Maslinic acid 
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Fig. 17(a): 7T9J C — Withaferin A 


Fig. 18(a): 7T9J_C — Withanone Fig. 16(b): 7T9J_C — Withanone 
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Fig. 19(a): 7T9J_C — Withanolide 
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Fig.5(a)-20(a): Crystal structure of the protein receptor-binding domains (7T9J_C) and ligand complexes of HCQ, p- 
Sitosterol, Betulinic acid, Desmodin, Foetidin, Glycyrrhetic acid, Kamolonol, Marrubiin, Maslinic acid, Oleanolic acid, 


Oleanane, Ursolic acid, Withaferin A, Withanone, Withanolide, and Zeylanone. 


Fig.5(b)-20(b): Hydrogen and hydrophobic interactions of all the compounds with the spike protein of the omicron variant 
(7T9J C). 


3.2 ADME studies & Drug likeliness prediction 


Besides, the calculating the binding energy and the 
molecular interactions between the protein and the ligand, 
compounds were further allowed for ADME analysis. 
These properties evaluate the ability of a compound to act 
as a drug, which follows Lipinski's Rule of Five. 


Lipinski's rule of five is a rule of thumb that describes the 
drug likeliness of a compound, with specific therapeutic 
properties that are likely to have chemical and physical 
properties that would make it an orally active drug. This 
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rule figures out a few of the molecular properties that are 
essential for the drug pharmacokinetics in the human body, 
such as Absorption, distribution, metabolism, and 
excretion. 


Lipinski's rule of five states that the chemical compound 
must follow the below-mentioned criteria to be considered 
a safe drug, 
1. Molecular mass > 500 Da 
2. A logarithm octanol-water partition coefficient 
(MlogP) > 4.15 
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3. Number of H-bond donors > 5 
4. Number of H-bond acceptors >10 


The violation of 2 or more conditions anticipates that the 
compound is not orally active [71]. Swiss ADME web 


International Journal of Environment, Agriculture and Biotechnology, 7(5)-2022 


server was used to study the ADME properties of the 
phytochemicals. Table 3 shows the ADME analysis data 
obtained from the SwissADME server for all the 15 
phytochemicals along with HCQ. 


Table 3: ADME properties of the few selected compounds 


Name of the Molecular MlogP H — Bond | H — | Violations Drug 
Phytochemical Weight («500 | (<4.15) | Acceptor | Bond Likeliness 
Da) (<10) Donor 
(<5) 

Hydroxychloroquine 335.87 2.35 3 2 0 Yes 
B- Sitosterol 414.71 6.73 1 1 1 Yes 
Betulinic acid 456.70 5.82 3 2 1 Yes 
Desmodin 382.41 2.06 6 1 0 Yes 
Foetidin 382.49 3.83 4 1 0 Yes 
Glycyrrhetic 470.68 4.87 4 2 1 Yes 
acid 

Kamolonol 398.49 3.00 5 1 Yes 
Marrubiin 332.43 2.76 4 1 0 Yes 
Maslinic acid 472.70 4.97 4 3 1 Yes 
Oleanolic acid 456.7 5.82 3 2 1 Yes 
Oleanane 412.73 9.05 0 0 1 Yes 
Ursolic acid 456.70 5.82 3 2 1 Yes 
Withaferin A 470.6 2.75 6 2 0 Yes 
Withanone 470.6 2.75 6 2 0 Yes 
Withanolide 470.6 2.75 6 2 0 Yes 
Zeylanone 374.34 0.69 6 2 0 Yes 


From Table 3, we can observe that all the 15 
phytochemicals followed the Lipinski rule of five. Out of 
which, 8 phytochemicals including Desmodin, Foetidin, 
Kamolonol, Marrubiin, Withaferin A, Withanone, 
Withanolide, and Zeylanone showed no violations. And 
the remaining 7 phytochemicals including B- Sitosterol, 
Betulinic acid, Glycyrrhetic acid, Maslinic acid, Oleanolic 
acid, Oleanane, and Ursolic acid showed one violation 
(MlogP < 4.15) in Lipinski rule of five. 


The data obtained from Table 3, shows that all the selected 
phytochemicals followed the Lipinski rule of five with not 
more than 1 violation, making them a suitable drug 
candidate for the treatment against the omicron variant. 


The current study was focused on identifying the 
inhibitory potential of the 15 selected phytochemicals by 
targeting the active sites of SARS CoV-2 omicron variant 
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spike protein. This protein was chosen as the druggable 
target site, since it plays a crucial role in the recognition, 
attachment, and fusion of the virus into the host cell. There 
are several chemotherapeutic drugs such as Remdesivir, 
Azithromycin, Amantadine, and many more that are being 
used to treat the infection. But, few of the drugs showed 
certain side effects like low blood pressure, vomiting, 
nausea, orthostatic hypotension, allergic reactions, and 
others. In addition to this, higher uptake of 
chemotherapeutic drugs might cause severe toxic effects 
on the human body that can lead to organ dysfunction. 


Hence, there is an urgent need for developing an effective 
drug or inhibitor against the omicron variant with the least 
toxic effects. As a consequence of that, these 
phytochemicals can be the best solution for the treatment, 
whose binding energies were found to be 1.5-1.8 times 
higher than the std. HCQ, following the Lipinski rule of 
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five criteria and are known to be safer with no or fewer 
side effects, thus intimating that these phytochemicals 
have the potential to form an antiviral drug or inhibitor 
against the SARS CoV-2 omicron variant. 


IV. CONCLUSION 


The mutating nature of the SARS CoV-2 has resulted in 
the emergence of omicron; another new variant with a 
higher fatality rate has created urgency for the 
development of new antiviral drugs or vaccines against the 
infection. In this regard, the research study was carried out, 
to identify the therapeutic potential of the phytochemicals 
from different medicinal plants against the spike protein of 
the omicron variant using Insilico approaches. This study 
showed that the binding energy of all the 15 
phytochemicals was higher compared to the std. HCQ with 
the greater potential of inhibiting the binding of the viral 
spike protein to the host cell receptor. In addition to this, 
the phytochemicals also fulfilled all the criteria of the drug 
likeliness and ADME studies. Therefore, the results 
obtained from the current study suggest that all these 15 
phytochemicals can be used as lead molecules for 
designing an effective drug against the omicron variant. 


However, the data obtained from the current study is based 
on the Insilico approach, additional studies have to be 
carried out with in vitro and in vivo conditions using 
animal models to check the feasibility of the compound. 
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